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ABSTRACT 



< 



. , Context. Barium is a key element in constraining the evolution of the (not well understood) r-process in the first galactic stars and 

\^ . currently the Ba abundances in these very metal-poor stars were mostly measured under the Local Thermodynamical Equilibrium 

(UTE) assumption, which may lead in general to an underestimation of Ba. 

Aims. We present here determinations of the barium abundance taking into account the non-LTE (NLTE) effects in a sample of 
extremely metal-poor stars (EMP stars): 6 tumoff stars and 35 giants. 
^ ■ Methods. The NLTE profiles of the three unblended Ban lines (4554 A, 5853 A, 6496 A) have been computed. The computations 

1 ^-^ ' were made with a modified version of the MULTI code, applied to an atomic model of the Ba atom with 31 levels of Bai, 101 levels 

lyj ' of Ban, and compared to the observations. 

d . Results. The ratios of the NLTE abundances of barium relative to Fe are slightly shifted towards the solar ratio. In the plot of 

[Ba/Fe] versus [Fe/H], the slope of the regression line is slightly reduced as is the scatter. In the interval -3.3 < [Fe/H] < -2.6, 
[Ba/Fe] decreases with a slope of about 1.4 and a scatter close to 0.44. For [Fe/H] < -3.3 the number of stars is not sufficient to 
J decide whether [Ba/Fe] keeps decreasing (and then CD-38:245 should be considered as a peculiar "barium-rich star") or if a plateau 

it" , is reached as soon as [Ba/Fe] x -1. In both cases the scatter remains quite large, larger than what can be accounted for by the 

^— ^ . measurement and determination errors, suggesting the influence of a complex process of Ba production, and/or inefficient mixing in 

*0 ' the early Galaxy. 
'NT ■ 

• Key words. Galaxy: abundances - Galaxy: halo - Galaxy: evolution - Stars: abundances - Stars: Mixing - Stars: Supernovae 

(N 
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^\ 1. Introduction stars obtained in the same conditions. On the other hand, from 

T . . ^ , T-o^ ,„ „ ^ similar LTE analysis, e.g. Burris et al. ( |2000| l, Honda et al. 

^ In the framework of the ESO large programme First stars , ^ai et al. mm found abundance ratios that, in the 

, ^ ; very high quality spectra of extremely metal-poor (EMP) giants ^^^^ metallicity, show a similar behaviour This general 

, and turnoff stars were obtained with the high resolution spec- ^^^^^ ^j^^ ^^^^^^ ^^^^ ^^^^^^ measurements 
trographUVES fed with the VLT. In this sample, 33 giants and Barklemetal. (2005) in a very large sample of 373 metal-poor 
18 turnoff stars are normal metal-poor stars (not cai^on-rich). ^^^^^^ ^^^^-^ ^-^^ ^ ^ 

Among them 22 giants and 10 turnon stars have [Fe/H] < -3. ^ n u • u j ■ » i . u 

rr,, . r , , . • ^ , Generally, the barium abundance m metal-poor stars has 

The main parameters of these stars can be found in Cayrel et , . . j j »u tt-t- • u *u 

1 J7v7v7m„i-- r„~. 1 Ti ,c ■ 1 i^oo-T i ill- been investigated under the LTE approximation. However, the 

al. 0004 "First stars V ) and Bonifacio et al. 0007 First stars xttt^c u • ^ . • i . i 

, --, „ „ i ^,,TT„N T • ^ — — - ^AK- , NLTE corrections can be very signmcant, in particular at low 

\S u u 'T h f metallicity (see Asplund, gOOll Short & Hauschildt, » they 

stars VIII ) the abundance of the neutron capture elements has , . , \, . ^, ^ ^, , ~ — c ^i, 

, ■• , • , , jTT-, • T T^T- 1 could be at least partly responsible for the large scatter of the 

been studied in the sample of EMP giants. From an LTE analy- > j f u i . . n- t 

. . , , „ , , „ , • , , , abundance ratios of the heavy elements at low metallicity. In 

sis. It has been found that, generally speaking, the abundance ra- j.j. • lui. j j i 

. ' , , '? . ■' 5„,„ , order to determine reliable trends and scatters, and to properly 

tios of these elements relative to iron ([X/Fe]) are very scattered, , u - . t i • ■ c f . 

., . ^ ^ ■ ■ , ,, ■ trace the history of nucleosynthesis processes in the hrst stars. It 

with sometimes a factor of more than one hundred between the • u- ui j • ui » c ^^^^rr^^ i 

, „ IS highly desirable to perform NLTE analyses, 
abundance ratios of two stars with equal [Fe/H]. Moreover, on ° i j- j j i j ■ ■ r , 

„ , . J, , 1 ^ , •, rr- /m There are several studics devotcd to thc determination of the 

average, [Ba/Fe] increases from about -1.5 to -1-0.2 while [Fe/H] • i j • i ^- /,rvoo^ 

c -, £ . T ^ -ru • . . „ r rn n NLTE barium abundancc in mctal-poor stars. Gigas (1988) per 

increases from -3.6 to -2.6. The important scatter of [Ba/Fe] „ , xttt-t- j i- . i ■ i- • f 



/ J rn rr- /.ti 1 i n j i n vc • i foiTOcd NLTE modeling of the barium lincs in the spcctrum of 

(and [Sr/Fe]) vs. [Fe/H] has been conhrmed by Bonifacio et al. , , j ^ ■ j r ^r t «i i i • 

Jt^- i^tt^.x 1 1 1 J 1 , r the mildly metal dehcient dwarf star Vega. Later, Mashonkina 

02OO 8 First stars XII ) where the abundances in the sample of „ /ir,r,^^ j . i- i /ir,r,r,x • j 

^■^ — , J . u J ■ 1 f . & Bikmaev (1996), and Mashonkina et al. (1999) estimated 

giants are compared to the abundances in a sample of turnoff , ■ ■ r , , j i j • 

° ^ ^ abundance correction for the cooler and more metal dehcient 



* Based on observations obtained with the ESO Very Large Telescope dwarfs (down to [Fe/H] = -2.5). This work was continued by 
at Paranal Observatory (Large Programme "First Stars", ID 165.N- Mashonkina & Gehren (2000), Mashonkina et al. (2003) and 
0276; PL: R. Cayrel. more recently by Mashonkina et al. (|2008> . In this last paper the 



2 



Andrievsky: NLTE Barium abundances in EMP stars 



authors study the abundance of the heavy elements in four very 
metal-poor stars with [Fe/H] w -2.5: one giant, HD 122563, 
and three turnoff stars. Moreover Short & Hauschildt (2006), 
performed a theoretical analysis of the non-LTE effects on the 
Ba lines {A = 4554, 5853, 6141, 6496 A) for metallicities be- 
tween [Fe/H] = -1 and -5. They used NLTE models for these 
computations, but an interesting point is that they could show 
that the NLTE profiles based on NLTE models differ negligibly 
from the NLTE profiles computed from LTE models (as is gen- 
erally done). 

We present here the results of a direct application of an 
NLTE analysis to the sample of extremely metal-poor stars 
(turnoff stars and giants) previously analyzed in the "First stars" 
programme. In our previous papers (Andrievsky et al. 2007, 
2008) we reported the sodium and aluminium NLTE abundances 
in these same stars. 

This sample covers the region of metallicities from [Fe/H] x 
-2.5 to-4.0. Such metal-poor stars have not been investigated up 
to now with the aim of deriving their NLTE barium abundances. 

2. Observations and reduction. 

The sample of stars and the observational data are the same as 
dicussed in previous papers : First Stars V, VII, VIII and XII 
(Cayrel et al.. [20041 Bonifa cio et al.. l2007t Frangois et al.. l2007l 
and Bonifacio et al.. 120081 ). In brief, the observations were per- 
formed with the VLT-UT2 and UVES (Dekker et al., |2000] l with 
a resolving power R ^ 45000. The signal-to-noise ratio in the 
spectra of the giants is very high: ~ 130 per pixel (with an av- 
erage of 5 pixels per resolution element) but it is lower for the 
turnoff stars (generally fainter), in particular in the /14554A re- 
gion where the strongest Ba line is located, at the very end of the 
spectra obtained with the "blue" camera. 

The spectra were reduced using the UVES context (Ballester 
et al. |200D] i within MIDAS. 

Table 1 gives the NLTE barium abundances for our sam- 
ple. They could be determined for only six (out of 18) turnoff 
stars, because even their strongest lines of Ba are very weak. 
The e(Ba)LTE value, given in column 6 for comparison, has been 
borrowed from Francois et al (2007) for the giants and Bonifacio 
et al. (2008) for the turnoff stars. In the last column is given the 
total number of Ba lines used for the computation. If n=3, all the 
three lines 4554, 5853, and 6496A have been computed, if n=2 
only 4554 and 6496 A could be used and if n=l only 4554A. For 
the giants, an "m" indicates that the star has been found "mixed" 
in Spite et al. (120051 ). In this table, for the determination of the 
relative LTE and NLTE values we have adopted for the solar 
values loge(Ba) - 2.17 (Asplund et al. 120051 ) at variance with 
Fran9ois et al. 120071 who adopted loge(Ba) - 2.13 (Grevesse & 
Sauval 120001). 

3. NLTE calculations. 

The NLTE profiles of the barium lines were computed using a 
modified version of the MULTI code (Carlsson 1986i ). The mod- 
ifications are described in Korotin et al. ( 1999). 

For these computations we used Kurucz's 1992 models 
(ATLAS9), after checking on some typical stars that the use 
of MA RCS models (Gustafsson et al. 120081) as in Cayrel et al. 
( 1200 4^ would not introduce any significant difference. 

As a rule, in the turnoff stars with [Fe/H] < -2, the barium 
abundance can be determined only from the resonance line at 
4554 A, but even this line becomes too weak at [Fe/H] < -3.0. 



In the giants, three lines can be traced down to [Fe/H] ^ -3.0, 
the resonance line and the two subordinate lines at 5853 and 
6496 A.; at still lower metalhcity only two lines (4554 and 6496 
A ) are still detectable. 



3. 1 . Atomic model of Ba 

Our barium model contains 31 levels of Bai, 101 levels of Ban 
with n < 50, and the ground level of Ba iii ion. In the detailed 
consideration we included 9 1 bound-bound transitions between 
the first 28 levels of Ban with n < 12 and I < 5. The other levels 
were used for the particle number conservation. The fine struc- 
ture of the levels 5d^D and 6p^P'' was taken into account; the 
other levels were treated as single. The corresponding Grotrian 
diagram is shown in FiglT] Only those transitions that were con- 
sidered in detail, are indicated. The energies of the levels are 
from Curry ( 120041 ). The oscillator strengths of the bound-bound 
transitions are from Wiese & Martin( 1980), Warner ( 119681 ). and 
Miles & Wiese ( 119691 ). For the transitions between the seven 
low-lying levels we used the data from Davidson et al. d 19921 ). 



90000 




Fig. 1. The Grotrian diagramme for Ban, only those transitions 
that were considered in detail, are indicated. 
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Fig. 2. Profile fitting for the barium lines in the solar spectrum 
withloge(Ba) = 2.17. 
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Table 1. Adopted model and NLTE barium abundance for our sample of stars. An asterisk after the name of the star means that the 
star is carbon-rich. In the last column is given the number of Ba lines used for the computation; the letter m indicates that the giant 
has been found "mixed" (see text). 





star 


TefF,(K) 


logg 


Vt,(kms ') 


[Fe/H] 


6(Ba)LTE 


e(Ba)NLTE 


[Ba/H]NLTE 


[Ba/Fe]NLTE 


Rem 


turnofF stars 




BS 17570-063 


6240 


4.8 


0.5 


-2.92 


-1.05 


-0.83 


-3.00 


-0.08 






CS 22948-093 


6360 


4.3 


1.2 


-3.30 


-0.90 


-0.83 


-3.00 


+0.30 






CS 22966-011 


6200 


4.8 


1.1 


-3.10 


-0.95 


-0.73 


-2.90 


+0.20 






CS 29506-007 


6270 


4.0 


1.7 


-2.91 


-0.55 


-0.23 


-2.40 


+0.51 






CS 29506-090 


6300 


4.3 


1.4 


-2.83 


-1.00 


-0.63 


-2.80 


+0.03 






CS 30301-024 


6330 


4.0 


1.6 


-2.75 


-0.85 


-0.53 


-2.70 


+0.05 




giants 


1 


HD 2796 


4950 


1.5 


2.1 


-2.47 


-0.48 


-0.46 


-2.63 


-0.12 


3, m 


2 


HD 122563 


4600 


1.1 


2.0 


-2.82 


-1.59 


-1.45 


-3.62 


-0.80 


3, m 


3 


HD 186478 


4700 


1.3 


2.0 


-2.59 


-0.50 


-0.45 


-2.62 


-0.03 


3, m 


4 


BD +17:3248 


5250 


1.4 


1.5 


-2.07 


+0.75 


+0.48 


-1.69 


+0.38 


3, m 


5 


BD -18:5550 


4750 


1.4 


1.8 


-3.06 


-1.67 


-1.45 


-3.62 


-0.56 


3 


6 


CD -38:245 


4800 


1.5 


2.2 


-4.19 


-2.82 


-2.50 


-4.67 


-0.48 


1, m 


7 


BS 16467-062 


5200 


2.5 


1.6 


-3.77 


< -2.76 


< -2.35 


< -4.52 


< -0.75 




8 


BS 16477-003 


4900 


1.7 


1.8 


-3.36 


-1.68 


-1.40 


-3.57 


-0.21 


2 


9 


BS 17569-049 


4700 


1.2 


1.9 


-2.88 


-0.55 


-0.50 


-2.67 


+0.21 


3, m 


10 


CS 22169-035 


4700 


1.2 


2.2 


-3.04 


-2.10 


-1.95 


-4.12 


-1.08 


2, m 


11 


CS 22172-002 


4800 


1.3 


2.2 


-3.86 


-2.90 


-2.60 


-4.77 


-0.91 


1, 


12 


CS 22186-025 


4900 


1.5 


2.0 


-3.00 


-0.85 


-0.72 


-2.89 


+0.11 


3, m 


13 


CS 22189-009 


4900 


1.7 


1.9 


-3.49 


-2.65 


-2.40 


-4.57 


-1.08 


1, 


14 


CS 22873-055 


4550 


0.7 


2.2 


-2.99 


-1.31 


-1.15 


-3.32 


-0.33 


3, m 


15 


CS 22873-166 


4550 


0.9 


2.1 


-2.97 


-1.54 


-1.32 


-3.49 


-0.52 


3, m 


16 


CS 22878-101 


4800 


1.3 


2.0 


-3.25 


-1.70 


-1.35 


-3.52 


-0.27 


2, m 


17 


CS 22885-096 


5050 


2.6 


1.8 


-3.78 


-2.75 


-2.50 


-4.67 


-0.89 


1 


18 


CS 22891-209 


4700 


1.0 


2.1 


-3.29 


-1.71 


-1.43 


-3.60 


-0.31 


3, m 


19 


CS 22892-052* 


4850 


1.6 


1.9 


-3.03 


+0.11 


-0.03 


-2.20 


+0.83 


3 


20 


CS 22896-154 


5250 


2.7 


1.2 


-2.69 


-0.05 


-0.13 


-2.30 


+0.39 


3 


21 


CS 22897-008 


4900 


1.7 


2.0 


-3.41 


-2.28 


-2.15 


-4.32 


-0.91 


1 


22 


CS 22948-066 


5100 


1.8 


2.0 


-3.14 


-1.95 


-1.75 


-3.92 


-0.78 


l,m 


23 


CS 22949-037* 


4900 


1.5 


1.8 


-3.97 


-2.42 


-2.25 


-4.42 


-0.45 


l,m 


24 


CS 22952-015 


4800 


1.3 


2.1 


-3.43 


-2.63 


-2.35 


-4.52 


-1.09 


1, m 


25 


CS 22953-003 


5100 


2.3 


1.7 


-2.84 


-0.22 


-0.23 


-2.40 


+0.44 


3 


26 


CS 22956-050 


4900 


1.7 


1.8 


-3.33 


-1.98 


-1.73 


-3.90 


-0.57 


1 


27 


CS 22966-057 


5300 


2.2 


1.4 


-2.62 


-0.73 


-0.73 


-2.90 


-0.28 


2 


28 


CS 22968-014 


4850 


1.7 


1.9 


-3.56 


< -3.20 


< -2.80 


< -4.97 


< -1.41 




29 


CS 29491-053 


4700 


1.3 


2.0 


-3.04 


-1.80 


-1.52 


-3.69 


-0.65 


3, m 


30 


CS 29495-041 


4800 


1.5 


1.8 


-2.82 


-1.34 


-1.06 


-3.23 


-0.41 


3 


31 


CS 29502-042 


5100 


2.5 


1.5 


-3.19 


-2.75 


-2.40 


-4.57 


-1.38 


1 


32 


CS 29516-024 


4650 


1.2 


1.7 


-3.06 


-1.83 


-1.45 


-3.62 


-0.56 


3 


33 


CS 29518-051 


5200 


2.6 


1.4 


-2.69 


-1.01 


-0.60 


-2.77 


-0.08 


3, m 


34 


CS 30325-094 


4950 


2.0 


1.5 


-3.30 


< -3.05 


< -2.85 


< -5.02 


< -1.72 




35 


CS 31082-001 


4825 


1.5 


1.8 


-2.91 


+0.39 


+0.05 


-2.12 


+0.79 


3 



3.2. The different constants 

A cause of uncertainty in the NLTE analysis of the barium spec- 
trum is the scarce information about the photoionization cross- 
sections for the different levels. For the majority of the levels we 
used the results of the scaled Thomas-Fermi method application 
(Hofsaess'1979'). For g-levels and levels with n = 11 we used the 
hydrogen-like approximation (Lang, 1998 , eq. 1.230). 

Effective excitation electron collisional strengths for the 
transitions between the first levels (6s^S, 5d^D and 6p^P*') were 
used following Schoening & Butler ( 1998). Experimental cross- 
sections for the transitions 6s^S - 7s^S and 6s-S - 6d^D were 
taken from Crandall et al. ( 1974). Collisional rates for the tran- 
sitions between sublevels 5d^D, 6p^P'^ and 7s^S, 6d^D, and be- 
tween 7s^S and 6d^D were estimated with the help of the cor- 



responding formula from Sobelman et al. (119811 ). For the rest 
of the allowed transitions, we used the van Regemorter ( l962| 
formula, while for the forbidden transitions the Allen's ( 1973] ) 
formula was used. 

The rate of the collisional ionization from the ground level 
of Ban was calculated using the corresponding formula in 
Sobelman et al. ( |1981| l. The use of experimental data of Peart 
et al. d 19891 1 and Feeney et al. (1972) give the same results. For 
the other levels, Drawin's ( 1961) formula was used. 

Inelastic collisions of barium atoms with hydrogen atoms 
may play a significant role in the cool star atmospheres. They 
have been taken into account through the formula of Steenbock 
& Holweger ( |1984| |. A correcting factor of 0.1 was derived as a 
result of the experimental fitting of the barium lines in solar spec- 
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-3 
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> 

m 
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0.25 



Fig. 3. Ratio S;/Bv (Source function over Planck function) vs. 
the optical depth for the models Tgff - 5250K, log g - 
1.5, [Fe/H] - -2.0 and -3.0. The dashed line represents the 
4554 A line, the solid line the 5853 A line and the dotted line 
the 6496 A line. 



trum. Collisions with atomic hydrogen for the forbidden transi- 
tions were not taken into account, but we checked that the result- 
ing uncertainty about the equivalent width of the lines, computed 
by the method of Takeda (1991), does not exceed 0.5%. 

The odd barium isotopes have hyper-fine splitting of their 
levels and thus several HFS components for each line. As was 
demonstrated by Mashonkina et al. (1999), a three-component 
structure is sufficient to describe the 4554 A barium line (see 
also Rutten 1978). This line was fitted in the solar spectrum by 
adopting the even-to-odd abundance ratio of 82:18 (Cameron 
I1982I I. Since metal deficient stars are supposed to include a 
larger fraction of the material synthesized in SNe II through 
the pure r-process, we used an even-to-odd ratio 50:50 for the 
stars. Nevertheless, we checked that the resulting difference in 
the equivalent widths of the 4554A line due to this change in 
the even-to-odd ratio is very small. (As a consequence, from 
our spectra it is not possible to estimate the abundance ratios of 
the different isotopes of Ba from our spectra, higher resolution 
would be necessary.) 

Radiative damping constants are from Mashonkina & 
Bikmaev ( 119961 ), Stark broadening parameters are from the 
VALDatabase ( |http ://ams . astro .univie .ac . at/~ vald ) . Since the 
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Fig. 4. NLTE corrections computed as a function of the temper- 
ature of the model, for several metallicities. The calculations 
are performed for log g - \5 and [Ba/Fe]=-0.5 (solid lines). 
However, the computations have been done also for [Fe/H]= - 
3.0 and [Ba/Fe] = 0.0, thus [Ba/H] = -3.0 (dotted Hnes). It can 
be seen that the correction is then very similar to the correction 
computed for [Fe/H] - -2.5 and [Ba/Fe] = -0.5 (thus also with 
[Ba/H]= -3.0). The main parameter that determines the NLTE 
correction, is indeed [Ba/H] and not [Fe/H]. 



classical van der Waals formula underestimates the effect of 
the interaction with neutral particles, the broadening parameters 
were found by comparing the observed lines in the solar spec- 
trum (Kurucz et al. 1984) with a synthetic spectrum computed 
with the solar model of Kurucz ( 1996 ). The best agreement was 
obtained for loge(Ba) - 2.17 (Fig. |2]i, in good agreement with 
the estimation of the solar abundance of barium of Asplund et 
al. (2005). 

The oscillator strengths of the lines components and the 
corresponding HFS shifts, as well as the broadening parame- 
ters are given in Table |2] For the Van der Waals constant C(, 
we have adopted the definition of Unsold ( 119681 eq. 82,47): 
Ce = Aw rV27r. 
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Fig. 5. NLTE corrections computed as a function of the temper- 
ature of the dwarf models, for several metallicities. All calcula- 
tions are performed with [Ba/Fe]=0.0 and log g - 4.5. 
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Fig. 6. LTE (dashed lines) and NLTE (solid lines) profiles calcu- 
lated for Teff = 4800 K, [Ba/Fe]=0.0 and log ^ = 1.5. 



Table 2. Parameters of the barium lines. C(, - Aw I2n, follow- 
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3.3. Non-LTE effects 

The profiles of the barium lines are rather sensitive to the NLTE 
effects, in particular in metal-poor stars with rather high temper- 
ature. The NLTE correction sign changes when the temperature 
and the metallicity change. Mashonkina et al. ( 1999) discuss in 
detail this behaviour in detail for turnofF metal-poor stars. Since 
the abundance of barium at a given metallicity is very scattered, 
it is important to note that the NLTE corrections are very sensi- 
tive not only to the metallicity of the model but also to the barium 
abundance itself (the equivalent width of the line). 

- The resonance line or the low excitation subordinate lines 
that are used for the determination of the barium abundance can 



show substantial nonequilibrium excitation effects. In a metal- 
poor atmosphere, since there are few electrons and thus few col- 
lisions, the radiative processes cause the departure from equi- 
librium in the atomic level populations even in the deep atmo- 
spheric layers. If the mean intensity of the radiative field in the 
line frequencies Jy exceeds the Planck intensity B^, then an en- 
hanced photoexcitation depopulates the lower atomic levels and 
overpopulates the upper ones (this "pumping" is described e.g. 
in Bruls et al. fT992l or Asplund [2005l l. The depth of the effective 
formation of the radiation in b-b transitions changes significantly 
with the metallicity (and the corresponding barium abundance) 
of the star (Fig. |3]l. This explains the dependence of the NLTE 
correction on [Fe/H] (Fig. |4] and |5]l. The lower the metallicity 
of the model, the deeper in the atmosphere is the line formation 
level (Fig.O. 

- Another mechanism causing the deviations from LTE is the 
so called UV overionization from the 6p level. In the atmosphere 
of turnofF or RGB stars, Ba II is the dominant species and gener- 
ally this implies that overionization does not significantly affect 
the line formation. But for extremely metal-poor stars, already at 
T"5()oo - 1 the mean intensity exceeds the Planck intensity at the 
frequency of the ionization threshold of the 6p level. As a result 
this level is significantly depopulated. 

These effects are very sensitive to the metallicity of the 
model, the physical conditions and the characteristics of the ra- 
diation in the layers where the barium lines are formed. In Fig. 
|3]the ratio of the source function to the Planck function is plot- 
ted for two gjfint stars with the same parameters Tgff - 5250 K, 
log § = 1.5, Vt = 2 km s"', [Ba/Fe] = -0.5, but with differ- 
ent metallicities: [Fe/H] - -2.0 and -3.0. The formation depth 
of the three lines used for the determination of the Ba abun- 
dance (A - 4554, 5853, 6496A) is indicated in the figure. It 
can be seen that Si /By is larger than 1 at the depth of the line 
formation for the three considered lines for [Fe/H]= -3.0. In 
this case, NLTE effects make the line weaker. The model with 
a higher metal abundance produces another picture. While two 
lines 5853 A and 6496 A are formed in that region where 
Si/By > 1, the resonance line is formed significantly higher in at- 
mosphere, and for this line S//Bv < 1. This line is made stronger 
as a result of the departure from LTE. 

As an illustration we show in Fig. |4] (giants), and|5](turnoff 
stars) how the NLTE correction depends on the metallicity for 
different effective temperatures. For the giant stars we used 
[Ba/Fe] = -0.5, while for dwarfs [Ba/Fe] = 0.0. These values 
are typical of the observations. (In turnoff stars if [Fe/H]= - 
3.0 and [Ba/Fe]= -0.5, the strongest barium line is very weak 
(W < 4mA) and since it is located at the end of the spectra it is 
generally not detectable). In fact the main parameter which de- 
termines the NLTE correction is [Ba/H], as can be seen in Fig.|4] 
Our corrections for the 4554 A line are in good agreement with 
those of Mashonkina & Bikmaev (1 19961 1 computed for dwarfs 
and [Fe/H] =-2. 

In giant stars (Fig.lUi, the NLTE correction for the two lines 
at 4554 and 5853 A can be positive or negative depending on 
the effective temperature. A similar conclusion was obtained by 
Short & Hauschildt (I2006I I. For the line 4554 A we found that 
LTE and NLTE profiles coincide at [Fe/H] =-2.5, while Short & 
Hauschildt (2006) found such a coincidence at [Fe/H]= -3.0. 

Our computations (NLTE line transfer with an LTE back- 
ground model) cannot be directly compared to those of Short 
& Hauschildt who used NLTE line transfer on an NLTE back- 
ground model. In spite of this, the predictions of the two are 
qualitatively similar, with moderate quantitative disagreement. 
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Part of the disagreement may, indeed, be due to the difference in 
the adopted background models, but part of it may also be due 
to the differences in the model atoms employed in the two com- 
putations. Some examples of the calculated profiles are given in 

Fig.ia 



3.4. Line profile fitting 

It should be stressed that the NLTE corrections given in Fig.|4]or 
|5] should not be used to determine a precise value of the NLTE 
abundance of barium in metal-poor stars, since this correction 
depends critically on several parameters ([Ba/H], Vt). The best 
way to determine the Ba abundance is to calculate the NLTE 
profiles directly, and to compare them to the observed profiles. 
Such an approach was used in this paper. The profile fitting for 
some stars is displayed in Fig.|7] 

The mean NLTE barium abundances of the stars are listed in 
Table [T] . The small difference between the [Ba/Fe] ratio found 
for HD 122563 in Mashonkina et al. ('2008') and in our paper, 
can be explained by the different values adopted for the grav- 
ity. When several Ba lines are detected in the spectra of a star, 
the abundances deduced from the different lines agree within 
+0. Idex. This value can be considered as the observational error. 
The total error is the quadratic sum of this error and of the error 
due to the stellar parameter uncertainties. The largest uncertainty 
arises from the uncertainties in the temperature and the gravity 
of the stars. We estimate that the total error on [Ba/Fe] is close 
to 0.2 dex. 




4553.5 4554.0 4554.5 

X(A) 



Fig. 7. Profile fitting for the barium lines in two giant stars 
(BS 17569-049 and CS 22891-209) and a turnoff' star (BS 17570- 
063). In this last case only the 4554A line is visible. 
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Ba (NLTE) 



-3.5 



-3 
[Fe/H] 



-2.5 



Fig. 8. [Ba/Fe] vs. [Fe/H] in our sample of stars computed with 
the LTE hypothesis (plot a) and taking into account the NLTE 
effects (plot b). The values of [Ba/Fe] have been computed with 
loge(Ba) - 2.17 (Asplund et al., 20051) at variance with Frangois 
et al. (I2007I I who adopted loge(Ba) - 2.13. The filled circles 
represent the turnoff stars, the open circles (in red in the elec- 
tronic version of the paper) the RGB mixed giants and the open 
squares (in blue) the unmixed RGB giants. The two ("pink") 
star-symbols represent the two carbon-rich stars we have mea- 
sured (see Table [TJ. The thick (green) lines represent the mean 
slopes computed between [Fe/H]= -3.6 and [Fe/H]= - 2.6. The 
slope of the regression line is slightly shallower when non-LTE 
effects are taken into account. In plot b) (NLTE) the thin (green) 
lines are drawn at a distance of two sigmas from the regression 
line. The two carbon-rich stars are outside the two sigmas limit. 
If we consider only the "normal" stars (not carbon rich) then 3 
stars are outside this limit on the "Ba-rich" side. 



4. Discussion 

4.1. Global trend 

The behavior of the ratio [Ba/Fe] vs. metallicity for the investi- 
gated stars is displayed in Fig. [8] 

The turnoff and the giant stars show a similar behaviour 
(as judged from the few dwarfs for which Ba was measurable). 
There is no difference in the behaviour of mixed giants (stars lo- 
cated above the Luminosity Function Bump : Spite et al., .20031 
120061) and the low RGB giants. Barium is not affected by mixing 
along the RGB. 
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In Fran9ois et al. ( 120071 1 it has been shown from an LTE 
analysis that in the interval -3.5 < [Fe/H] < -2.6, the 
abundance of barium decreases rapidly with the metallicity, 
in agreement with previous works in the literature. Our NLTE 
analysis of the barium lines in the same extremely metal-poor 
star spectra has led to some slight shifts of the derived individual 
LTE [Ba/Fe] ratios toward the solar ratio, but the main tendency 
is only slightly changed. 

Globally the results may be summarised as follows (see 
FigM- 

Below the metallicity [Fe/H] w -3.0, the mean value of 
[Ba/Fe] is negative as found in Frangois et al. ( 120071 ). 

The slope of the NLTE value of [Ba/Fe] vs. [Fe/H] is a little 
shallower than in LTE: the slope of the regression computed in 
the interval -3.6 < [Fe/H] < -2.6 is about 1.14 + 0.12 in NLTE 
rather than 1.40+0.18 under LTE assumption (Fig.[8]l. This small 
variation of the slope was expected by Mashonkina et al. (120081) . 
The slope is compatible with a secondary production of Ba. 

In the region -3.6 < [Fe/H] < -2.6, the scatter around the 
regression line of the [Ba/Fe] vs. [Fe/H] plot (Fig. [8]) is a lit- 
tle smaller (0.44) for the NLTE values than for the LTE values 
(0.52), but it remains very large, much larger than observed for 
the abundances of other elements (where it is generally smaller 
than 0.2 dex). This scatter strongly exceeds what is expected 
from the measurement and determination errors, and illustrates 
the decoupling between iron and Ba nucleosyntheses. This point 
remains valid when using Mg (rather than Fe) as a reference el- 
ement, suggesting also a decoupling between the production of 
Mg (hydrostatic burning) and Ba (r-process). 

In the extremely low metallicity domain ([Fe/H] < -3.5), 
the mean behaviour of [Ba/Fe] is not clear. At this very 
low metallicity, the barium abundance could be estimated in 
only three "normal" stars. (CD -38:245, CS 22172-002 and 
CS 22885-096). In these stars [Ba/Fe] is around -0.7 suggesting 
that a "plateau" of [Ba/Fe] is reached. But for two other stars, 
only limits could be measured and for one of them it is rather 
low: [Ba/Fe] < -1.37 . 

4.2. A complex trend ? 

At very low metallicity ( [Fe/H] < -2.5), the behaviour of 
[Ba/Fe] with [Fe/H] (Fig.|8]) can be interpreted in three different 
ways: 

1. Our Fig. [8] can be interpreted, as proposed by Frangois et 
al. (i2007 ), by a very large scatter of [Ba/Fe] in the region 
-3.2 < [Fe/H] < -2.8: a factor of about 1000 from star to 
star, and (for [Fe/H] < -3.3) a plateau close to the value 
[Ba/Fe] = -1. 

This interpretation was suggested by Fran9ois et al. ( I2007I I 
after merging their measurements with those of Honda et 
al. ( |2004t . But it seems that there is a systematic differ- 
ence between these two sets of measurements: in Fig. 10 of 
Franfois et al. the data points from Honda et al. seem to be 
systematically ofFsetted in [Fe/H] and [Ba/Fe] compared to 
the points of Frangois et al. (2007). The difference is mainly 
due to the adoption of a lower microturbulent velocity by 
Honda et al. (2004). This systematic shift could be respon- 
sible, at least partly, for the very large scatter without any 
defined trend, suggested by Franfois et al. ( 120071 ) in the re- 
gion -3.2 < [Fe/H] < -2.8. 



2. When only our measurements are taken into account (Fig.[8]l, 
and as a consequence, all the stars in the graph are studied in 
a very homogeneous way, it becomes clear that in the interval 
-3.6 < [Fe/H] < 2.5, [Ba/Fe] increases with [Fe/H] with a 
slope close to 1.14 (when the non-LTE effects are taken 
into account). 

A problem is to know what happens at lower metallicity. 
For [Fe/H] < -3.6 the decrease of [Ba/Fe] could continue 
and then 5 stars w ould seem Ba-rich with the definition 
[Ba/Fe] > [Ba/Fe] + 2cr : 

- the two carbon rich stars (CS 22892-052 (Sneden et al. 
[19961 and CS 22949-037 Depagne et al..2002 ), 

- the well known r-rich star CS 31082-001 (see Cayrel et 
al.,|200l] Hill etal., 1 20021 , 

- CD-38:245, the most metal-poor star of the sample. In 
this star Frangois et al. measured an upper limit of the 
europium abundance: [Eu/Fe] < +0.38. According to 
the definition (Barklem et al.. 120051 Beers and Christlieb, 
120051 ) of the r-I class of "moderately r-rich stars" (+0.3 < 
[Eu/Fe] < +1.0) the upper limit of the Eu abundance 
may bring CD-38:245 in this class. 

- a turnoff star, CS 22948-093. In this kind of star- it would 
be impossible to measure the europium abundance even 
if it is r-rich with [Eu/Fe] - 1 . 

All these stars are or could be peculiar (two of them are well 
known very r-rich stars). 

However if we suppose that in "normal" stars [Ba/Fe] de- 
creases continuously in the interval -4.2 < [Fe/H] < -2.5 
with a uniform scatter (sigma=0.44), then 4 stars (over 38) 
are outside the two sigma limits (3 on one side and 1 on 
the other side). This is more than what it is expected for a 
Gaussian population (only 5%). 

3. A third interpretation would combine a decrease of [Ba/Fe] 
only in the range -3.3 < [Fe/H] < -2.5 and then a plateau 
which would determine the mean value of the ratio [Ba/Fe] 
in the early Galaxy: [Ba/Fe] w -0.7 from this NLTE analy- 
sis. The value of this plateau would define the yields of bar- 
ium (relative to iron) in the massive primitive supernovae. If 
we consider only the stars with [Fe/H] < -3.7, the scatter 
around [Ba/Fe] = -0.7 is rather small, but there are 2 stars 
with upper limits that increase this scatter. If a plateau ex- 
ists, it could suggest a production of Ba independent of the 
metallicity in the primitive supernovae. 

The analysis of a larger sample of EMP stars with [Fe/H] < -3.3 
would obviously help in the interpretation of the results. 

5. Conclusion 

We present here a homogeneous determination of the abundance 
of Ba in a sample of extremely metal-poor stars (turnoff and gi- 
ant stars) taking into account the non-LTE effects. 

There is a good agreement between the abundances of turn 
off stars, unmixed giants (low RGB) and "mixed" giants (stars 
located in the HR diagram higher than the Luminosity function 
bump). In the giant stars the abundance of barium does not de- 
pend on the deep mixing inside the star. 

The NLTE abundances are free from the approximations 
made in the LTE determinations. The trend of [Ba/Fe] versus 
[Fe/H] is refined : less scatter, and a better defined and shallower 
slope. The slope is compatible with a secondary process. The 
general behaviour would be the same if Mg (rather than Fe) was 
the reference element. 
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The scatter of [Ba/Fe], although reduced, is large, and can- 
not be explained by the determination errors: it suggests the 
influence of several factors. The contribution of one (or more) 
additional "r" process has been invoked in the literature (see 
Francois et al. and references therein 120071 Travaglio et al., 
120041 Wasserburg & Qian, I2008I I. Also, an extremely efficient 
galactic mixing is excluded. 

The behaviour of [Ba/Fe] with [Fe/H] (Fig.O could be de- 
scribed in two different ways. For a decreasing metallicity : 

1. the mean value of [Ba/Fe] decreases linearly with a large 
(uniform) scatter and a few outliers 

2. the mean value of [Ba/Fe] decreases linearly in the interval 
-2.5 > [Fe/H] < -3.6, and then reaches a plateau. 

The number of observed stars at very low metallicity is ob- 
viously too low to reach any firm conclusion about a possi- 
ble plateau, that could suggest a production of barium at low 
metallicity independent of the metallicity. Surveys detecting ex- 
tremely metal-poor stars, and the spectroscopic observations of 
the detected stars, are required. 
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